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Quenched Hadron Spectroscopy with Improved Actions 

Wolfgang Bock^ 

'^Center for Computational Physics, University of Tsukuba, 305 Ibaraki, Japan 

A variety of different combinations of improved gluon and fermion actions is tested for quenched hadron 
spectroscopy. 



1. INTRODUCTION 

The improvement of the lattice QCD action has 
become recently a main focus of research in lattice 
field theory. Improved actions for lattice QCD are 
designed to remove lattice-spacing artifacts in nu- 
merical estimates of physical observables. This 
allows to conduct the numerical calculations on 
much coarser and smaller lattices and thereby 
to substantially reduce the simulation costs. In 
this contribution a variety of different combina- 
tions of gauge and fermion actions is tested for 
quenched hadron spectroscopy. For the gluon ac- 
tion we used the tadpole improved Liischer-Weisz 
(T-LW) |l|-|l and the Iwasaki (IW) action [|. 
The LW action has been obtained by extending 
Symanzik's perturbative improvement program 
to lattice gauge theories, demanding 0(a^) im- 



provement of on-mass-shell quantities m. The 
use of tad-pole improvement is essential to en- 
hance the convergence of the perturbative expan- 
sion and extend perturbation theory to larger dis- 
tances The IW action on the contrary is a 
renormalization group improved action. Scaling 
tests carried out recently for both models have 
been very encouraging The fermion ac- 

tions which we shall consider in this contribution 
are the Sheikholeslami-Wohlert (SW) action [Q, 
the D234 action Q and for comparison also the 
standard Wilson action (W). Symanzik improve- 
ment of the Wilson action leads to the SW action 
which eliminates the cut-off effects in on-mass- 
shell quantities to order 0{a). With the D234 
action one attempts to go even one step beyond 
the 0{a) improvement and to eliminate partially 
some of the cut-off effects to O(a^). Also in the 
case of the improved fermion actions it is essential 



to couple the fermions to a "tadpole improved" 
(T) lattice link variable We computed the 
hadron spectrum for the following five combina- 
tions of improved gluon and fermion actions: T- 
LW&W, T-LW&T-SW, T-LW&T-D234, IW&T- 
SW and IW&T-D234. 

2. DETAILS OF THE SIMULATION 

In order to be able to compare the results ob- 
tained with the different gluon actions we have 
conducted the simulations for the T-LW and IW 
actions at the critical (3 of the deconfining phase 
transition for A^t = 2 ^ . We performed the simu- 
lations on a 6^16 lattice using a standard Hybrid 
Monte Carlo algorithm. For the determination 
of the hadron spectrum we generated about 500 
configurations, which were separated by 50 tra- 
jectories. For the inversion of the fermion ma- 
trix we used the BICG75 algorithm which differs 
from the CG algorithm by the insertion of a 75 
in the inner products. We find that the BICG75 
is, within the interval = 0.7-1.0, by a fac- 
tor 3-4 times faster than the BICGstab and by 
a factor 6-8 faster than the CG algorithm. The 
BICG75 algorithm is not guaranteed to converge 
as a division by zero can occur during the iter- 
ation process. We have however never encoun- 
tered a breakdown of the algorithm although we 
performed in total about 2.5 10^ inversions of the 
fermion matrix. 

For the determination of the hadron spectrum 
we implemented the standard correlation func- 
tions for the TT, p, flo, ai and 61 mesons, the nu- 
cleon (N) and the A, using local operators for the 
quark field, both at the source and the sink. The 
meson and baryon propagators have been fitted 
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Figure 1. The mass ratios mN/rnp (a) and 
m\/mp (b) as a function of m-^lmp. The T- 
LW&W, T-LW&T-SW, T-LW&T-D234, IW&T- 
SW and IW&T-D234 data are marked by the 
filled triangles, open circles, crosses, open trian- 
gles and squares. The horizontal error bars are 
dropped as they are smaller than the symbol size. 

with AM[exp(— mAfi) + exp(— mM(T — t))] and 
j4b exp(— mst). We varied tmin systematically 
and averaged the masses when it was difficult to 
pick out the best among equally good fits. The 
mass spectrum has been determined in the in- 
terval = 0.7-1.5, where finite size effects are 
very small. All correlation functions gave a very 
clear and stable signal in this pion mass range, 
expect for the ag, a\ and h\ channels. We dis- 
covered however that the signal in those chan- 
nels improves significantly at larger pion masses, 
and therefore extended the calculation of the ao, 




0.5 1 1.5 2 

Figure 2. The masses rrip, m^r and mA as a 
function of m\ for the IW&D234 action. 

a\ and h\ correlation functions with the IW&T- 
D234 action, and for comparison also with the 
IW&W action, to somewhat larger pion masses, 
TOtt — 1.5-2.0. Apart from the hadron mass spec- 
trum we also computed the string tension a which 
we determined from the exponential fall-off of the 
Polyakov loop correlation function. 

3. RESULTS 

In fig. 1 we have plotted the mass ratios 
mj^/nip (a) and nii^/mp (b) as a function of 
Tfi^jmp. The filled circles to the left correspond 
to the physical point and the filled squares to the 
right to the heavy quark limit, where mT^/mp = 1 
and m-^lrrip — m/\/mp — 3/2. The solid lines 
connecting the two points represent the result of 
a phenomenological quark model formula derived 
in ref. Q . The graphs show that the points which 
correspond to the combinations T-LW&T-SW 
(open circle), T-LW&T-D234 (crosses), IW&T- 
SW (open triangles) and IW&T-D234 (squares) 
cluster in a narrow strip which is located in both 
cases above the solid line. The large gap be- 
tween the standard Wilson action data (filled tri- 
angles) and all other data reveals very clearly the 
dramatic effect of the fermion improvement. A 
naive linear extrapolation of the T-LW&T-SW, 
T-LW&T-D234 and IW&T-D234 data in fig. la 
down to the physical value of m^/mp leads to a 
totv /rup value which is not too far off from the ex- 
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Figure 3. Mass ratios for the various actions. The 
symbols have the same meaning as in fig. 1. 

perimental point (fiUed circle) . The extrapolation 
of the IW&T-SW data on the other hand results 
in a value that is significantly larger indicating 
that the cut-off effects are bigger for this partic- 
ular action. We checked that our data obtained 
with the improved fermion actions are, within er- 
ror bars, consistent with the standard Wilson ac- 
tion results which have been obtained during the 
past years on a much finer and bigger lattices. 

Using chiral perturbation theory we extrapo- 
lated the hadron masses m^, ma„, nia^, rrih-^, rriN 
and toa to the chiral limit, TOtt — > 0. To this 
extend we fitted the meson mass data rriM with 
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(0) 
M 



C2ml, 



where m^^j designates the mass in 



the chiral limit and M — p, qq, ai,bi in our case. 
For the nucleon it is essential to include also a 
term in the chiral fit which is cubic in the pion 
mass, m^'' +C2m'^ -j-cam^. Our results show that 
the coefficient C3 depends strongly on the particu- 
lar action. It is largest for the three combinations 
T-LW&T-SW, T-LW&T-D234 and IW&T-D234 
and roughly by a factor two smaller for the other 
In the case of the A we fitted 
as a cubic fit 



two actions 
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the mass data with m 
leads to a C3 value which is consistent with zero 
|lO[ . As an example, we plotted in fig. 2 the mp, 
rriN and uia data for the case of the IW&T-D234 



Table 1 


IW&W 


IW& 
T-D234 


experiment 


(0) / (0) 
rriai /mp 

(0) / (0) 
(0) / (0) 


2.37(15) 
2.37(17) 
2.28(14) 


1.81(8) 
1.70(10) 
1.90(10) 


1.600(52) 
1.279(2) 
1.601(13) 



are represented in this graph by the solid lines. 
The ratios of the extrapolated masses, 
and rri'^ /mf^ are displayed in fig. 3 
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along with ^Ja/m}p' . The horizontal lines rep- 
resent the experimental values. We find that 
the T-LW&T-SW (open circle), T-LW&T-D234 
(crosses) and IW&T-D234 (squares) mass ratios 
are very close to the experimental values. The 
gap is slightly bigger for the IW&T-SW data 
(open triangles), and largest for the standard 
fermion action data (filled triangles). The re- 
sults for the mass ratios niao lrn"p\ m^ai I'rn'p^ 



and m'^^ /mf^ 



action as a function of i 



The results of the fits 



are summarized in table 1 which 
also shows that the mass ratios obtained with the 
improved fermion action are much closer to the 
experimental values than the ones obtained with 
the standard fermion action. 
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